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ABSTRACT 

Afterglow observations are commonly used to determine the parameters of GRB explosions, the 
energy E, surrounding density n, post-shock magnetic field equipartition fraction es and electron 
equipartition fraction e e , under the frequently made assumption that the efficiency of electron "injec- 
tion" into relativistic shock acceleration is high, i.e. that the fraction / of electrons which undergo 
acceleration is / ~ 1. We show that the value of / can not be determined by current observations, since 
currently testable model predictions for a parameter choice {£" = E/f, n' — n/ f, e' B = feB, e' e = /e e } 
are independent of the value of / for m e /m p < / < 1. Current observations imply that the efficiency 
/ is similar for highly relativistic and for sub relativistic shocks, and plausibly suggest that / <~ 1, 
quite unlike the situation in the Crab Nebula. However, values m e /m p < / < 1 can not be ruled 
out, implying a factor m e /m p uncertainty in determination of model parameters. We show that early, 
< 10 hr, radio afterglow observations, which will be far more accessible in the SWIFT era, may provide 
constraints on /. Such observations will therefore provide a powerful diagnostic of GRB explosions 
and of the physics of particle acceleration in collisionless shocks. 



Subject headings: acceleration of particles — gamma-rays: bursts and theory 
— shock waves 



synchrotron radiation 



1. INTRODUCTION 

Synchrotron emission by shock accelerated particles is 
central to our understanding of explosive, high energy 
astrophysical phenomena, such as supernova remnants, 
jets from AGN and quasars, plerionic nebulae, and 7- 
ray burst (GRB) afterglows. GRB afterglows have pro- 
vided an unprecedented opportunity for diagnosing the 
blast wave and attendant shock acceleration, because 
their brevity in the observer's time frame and ultrahigh 
Lorentz factors allow rapid evolution of the synchrotron 
spectrum which can be observed over a wide span of 
wavelength regimes in real time. 

The afterglow radia tion of 7-ray bursts (GRBs) (e.g. 
iKulkarni et al. 1 |2000|1 is naturally explained as due to 
synchrotron emission of electrons accelerated in relativis- 
tic collisionless shocks driven by the GRB explosion into 
the medi um surrounding the GRB progenitor (for re - 
views see lPiran ll200rl iMeszaros 112002 IWaxman 1 120031. 
The energy released in the explosion leads to the forma- 
tion of a diverging shock wave, which propagates into 
the ambient plasma. At sufficiently late time (at times 
much longer than the burst duration) all the explosion 
energy is carried by the shocked ambient plasma (vanish- 
ingly small fraction of the energy remains in the ejecta 
produced by the explosion). The radiation is believed to 
be produced by electrons of the ambient medium, which 
are accelerated to high energy as they pass through the 
diverging shock. 

The dynamics of a spherical shock wave is determined 
by the explosion energy E, and by the surrounding 
medium number density n. If the initial GRB outflow 
is jet-like, an additional parameter, the jet opening an- 
gle 9j, is required in order to specify the flow. For a given 
shock dynamics, the luminosity and spectrum of emitted 

Physics Department, Ben-Gurion University, Beer-Sheva 
84105, Israel; eichler@bgumail.bgu.ac.il 

2 Physics Faculty, Weizmann Institute, Rehovot 76100, Israel; 



radiation are then determined by the fractions es and 
e e of shock thermal energy carried, respectively, by mag- 
netic field and electrons, and by the shape of the electron 
distribution function. The fraction of explosion energy E 
converted to thermal energy in the shock is determined 
by the hydrodynamics and is of order unity. The electron 
and magnetic field energy densities are therefore propor- 
tional to e e E and esE respectively. The electron distri- 
bution function is commonly assumed to be a power law 
of index p = — dlnn e /dlne e , where e e is the electron 
energy, above some minimum energy s e o (We use e to 
denote single particle energy, and E to denote the total 
flow energy). 

The processes of magnetic field generation and electron 
injection in collisionless shocks are not understood from 
basic principles and cb, e e and e e o cannot at present be 
determined theoretically. Rather, they are treated as free 
parameters of the model, constrained by observations. It 
is important to note here that e e o is, in general, an inde- 
pendent parameter of the model. It is a function not only 
of e e , but also of the fraction / of electrons assumed to 
be accelerated to beyond e e Q. It is commonly assumed, 
however, that / = 1, in which case e e o is uniquely deter- 
mined by the other model parameters. 

The observed afterglow synchrotron spectra constrain 
p, as is well known, and the break in t he power law de cay 
of afterglow flux is widely believed ( Rhoads 1 11999ft to 
establish the opening angle 6j in terms of E/n. Under the 
assumption that all the ambient electrons were injected 
to beyond e e o, i.e. / = 1, the remaining four parameters, 
{es, e e , n, E}, are fixed by the four observables v m (the 
frequency of maximum intensity), F m (the intensity at 
v m )i v coo i (the synchrotron cooling break), and v a (the 
self-absorption frequency). 

The results that have been deduced about p, cb, and 
e e , are consistent both with current knowledge and cur- 
rent ignorance about shock acceleratio n: In bursts where 
p can be determined accurately (e.e. IWaxman I Tl997a: 
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[Galama et al Ill998at iFrail. Waxman fc Kulkarni 11200(1 
Istanek et al. 1119991) v = 2.2 ± 0.1 is inferred. This value 
is consistent with the theoretical value of p derived for 
test particle acceleration in relativistic shocks via the 
first order Fermi mechanism, assuming isotropic diffu- 
sion of particles in momentum sp ace, y = 2.22 ± 0.02 ob- 
tained in numerical calculations dBednarz fc Ost rowski I 
Il998t iKirk et al. 1 120001 lAchterberg et al. 1 12001|) . and 
p = 20/9 obtained by a m ore recent analytic analysis 
l|Keshet fc Waxman 11200^1 . This value of p is not con- 
sistent with test particle results for large angle scatter- 
ing in relativistic shocks, which produce very hard spec- 
tra. It is, however, consistent with the value expected 
in the 100 Mev -10 GeV ran ge by non-linear theory 
for cosmic ray-mediated shock (Ellis on fc Eichler"lll985t 
Ell ison & Double 2002). Despite the agreement of the 
observed and theoretically derived values of p, assum- 
ing isotropic diffusion, it should be kept in mind that 
questions remain about diffusive shock acceleration, par- 
ticularly in regard to relativistic generalization and in 
regard to electron injection, and that there are alterna- 
tive acceleration processes (e.g. [Axons fc TavanillT994t 
INishikawaet al. H2004HHededal Il200^ ~ 

It is natural to hope that the values of eg, e e 
are universal since they are determined by the micro- 
physics of the collisionless shock. The constancy of 
p and of e e among different bursts is strongly sup- 
ported by observations. Universal values of p and e e , 
p k, 2 and e e ~ 0.1, typically inferred from most op- 
tical afterglows, are also inferred from the clustering 
of explosion energies iFrail et al. 1 1200 1|) and from X- 
ray afterglow luminosity 3 (iFreedman fc Waxman 1 1200 It 
IBereer . Kulka rm fc Frail II2003|) . The value of 6r is less 
well constrained by observations. However, in cases 
where es can be reliably constrained by multi waveband 
spectra, values close to equipar tition are inferred (e.g. 
IFrail. Waxman fc Kulkarni 112000^) . Such high values for 
tB and e e are remarkable and beg for an explanation. 
The magnetic field required for allowing electron accel- 
eration and emission of synchrotron radiation may con- 
ceivably be produced in the collisionless shock driven by 
the GRB explosion by Weibel instabilities or the lik e (e.g. 
Blandford fc Eichler Ill987t IGruzinov fc Waxman"1ll999t 
Medvedev fc Loeb 111999ft . or it may be that the acceler- 
ated particles mix with the magnetic field of the fireball 
itself. 

No l ess surprising is the conclusion by W axman"! 
(1997b), that e e o is close to jnipC 2 and that the low fre- 
quency radio spectra imply that there are relatively few 
electrons in the decade or two just below e e o- Had the 
electrons been picked up by shock acceleration at some 
much lower energy than jm p c 2 , the power law spectrum 
imparted by the shock acceleration would have extended 
down to much lower energies, and only a small minor- 
ity of them would have made it to 7m p c 2 or higher. In 
the case of the Crab nebula, for example, which contains 
perhaps the best-studied relativistic shock wave, this is 
indeed the case: most of the electrons in the nebula emit 
in the radio, and probably have Lorentz factors of order 

3 Apparently deviant v alues of p JChevalier fc Lil 119991 : 
IPanaitescu fc Kumar! 120(13) are inferred based on light curves, 
rather than spectra, and are sensitive to model assumptions (e.g. 
they depend on the assumed radial dependence of the ambient 
medium density). 



10 2 , which is many orders of magnitude lower than 'ym p c 2 
and even about a factor of 10 2 below 7m e c 2 . More will 
be said about his below. While this paper is not aimed at 
explaining this gaping difference between afterglows and 
the Crab Nebula, it motivates us to check the assumption 
that / = 1 in the case of the former. 

In any case, we are unable to determine from basic 
principles the efficiency of electron "injection" to beyond 
some threshold energy well beyond jm e c 2 . Even when 
the number of electrons beyond some injection threshold 
e e o is known, we are unable to determine theoretically 
the fraction / of total electrons that these high energy 
electrons represent. It is conceivable that a large fraction, 
1 — / ~ 1, of the electron population do not participate 
at the acceleration process and remain well below e e Q. 
This is discussed in § In § we discuss observational 
signatures of the existence of such non-injected thermal 
electrons in GRB-induced blast waves. Our main re- 
sults and their implications are summarized in § We 
discuss both the implications to GRB phenomenology 
and the implications for the theory of collisionless shock 
acceleration, in particular in the context of constraints 
imposed by observations on astrophysical systems other 
than GRBs. 

2. MODEL PREDICTION DEGENERACY 
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Fig. 1. — A schematic representation of the post shock electron 
distribution, for a relativistic shock of Lorentz factor 7 (or sub- 
relativistic shock of velocity v). Scattering of electrons streaming 
toward the shock with Lorentz factor 7 (or velocity v) results in 
post shock "thermal" energy of ~ 7m e c 2 (or ~ m e v 2 ). A fraction 
/ of the electrons is assumed to be injected to the acceleration 
process, which significantly increases the average energy of these 
electrons, to ~ e e o, and produces a power-law distribution at e e > 
e e o- As we show here, afterglow observations imply e e o ~ ■ym p c 2 
in the relativistic phase and e e o ~ m p v 2 in the sub-relativistic 
phase, but do not allow to determine /. Afterglow observations 
also require an electron number density that increases with energy 
sufficiently fast, q = rflnn e /rfrne e > 0, over ~ 1.5 decades of 
energy below e e0 IWaxman I l ill997bl) . 

To clarify the issues involved in the electron "injection" 
problem let us consider the situation illustrated in fig- 
ure n which may arise for a relativistic shock propagat- 
ing with Lorentz factor 7^1 (or sub-relativistic shock 
propagating with velocity v <C c) into a cold plasma of 
protons and electrons (as may be the case for a shock 
driven by a GRB explosion into the ISM). In the shock 
frame, a cold stream of protons and electrons approaches 
the shock with Lorentz factor 7 (velocity v). The parti- 
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cles are being scattered at the shock front, resulting in a 
velocity distribution which is close to isotropic behind the 
shock, thus converting a large fraction of the kinetic en- 
ergy of the incoming flow to thermal energy. Isotropiza- 
tion of the electron and proton incoming flow would lead 
to a post shock proton "temperature" T p ~ jm p c 2 (or 
T p ~ nipV 2 ), and to a post shock electron "temperature" 
T e ~ 7TO e c 2 <C T p (or T e ~ m e v 2 ). In order for the 
electrons to gain a significant fraction of the post-shock 
thermal energy, some process must couple them to the 
protons, and accelerate them to energy ^> T e . This pro- 
cess is yet unknown, and we can not determine based on 
theoretical considerations what fraction of the electrons 
are being accelerated. Thus, in addition to e e , the ac- 
celeration process must be described by (at least) one 
additional parameter, the fraction / of accelerated elec- 
trons. We show here that afterglow observations imply 
£e0 ~ jm p c 2 in the relativistic phase and e e o ~ m p v 2 in 
the sub-relativistic phase, but do not allow one to deter- 
mine /. 

A s pointed out in iWaxman I l)1997b[K the energy dis- 
tribution of electrons below the characteristic accelera- 
tion energy, e e o, is constrained by radio observations. 
The slope of the radio afterglo w spectrum observed i n 
several cases (e.g. fig. 1 of iGalama et aT~lll998aj) . 
f u oc ;W 3 , is consistent with that expected for radia- 
tion emitted by electrons at e e o at frequencies well below 
their characteristic synchrotron frequency, which is some- 
what below optical at the observation time (typically of 
order days). In order for the emission from lower en- 
ergy electrons, at e e < e e o, not to modify this spectrum, 
q = d\nn e / d\ne e > —1/3 is required. A somewhat more 
stringent constraint may be obtained from the require- 
ment that the self-absorption optical depth produced by 
these electrons not be large enough to affect the observed 
self-absorption frequency i/ a , q = dlnn e /dlne e > 2/3, 
(provided of course that v a is unambiguously established 
by electrons at e e o). The uncertainty in the value of v a 
(and the values of other characteristic frequencies) de- 
termined by observation s, relaxes the latt er constraint to 
q = d\nn e /d\ne e > (IWaxman 111 997b|) . These state- 
ments hold for about 1.5 decades of energy below e e0 , 
i.e. the energy range over which the electron distribution 
would affect the radio emission that has been observed 
to date. Determining the electron spectrum below e e o is 
of interest for both GRB phenomenology and particle ac- 
celeration theory, suggesting that more careful analysis 
of radio spectra are warranted for obtaining better con- 
straints on q. This is, however, beyond the scope of the 
current paper, where we focus on the injection efficiency, 
/■ 

In what follows we discuss the degeneracy of afterglow 
model predictions, showing that the predictions obtained 
assuming / = 1 for some choice of model parameters, 
{E, n, es, £e}, are the same as those obtained for any 
value of /, m e /m p < f < 1, and {£" = E/f,n' = 
n/f, e' B = fcB,£' e — /e e }- In S 12.11 we discuss the hydro- 
dynamics of the flow, and in § 12. 21 we discuss emission of 
radiation. 

2.1. Hydrodynamics 

The apparent physical size of the radiation emitting 
region has been determined in several cases. During the 



relativistic stage of shock expansion, the size of the emit- 
ting region has been determined directly through the 
observati on of the suppres s ion of diffractive radio scin- 
tillat ion l|Goodman I Il997t IWaxman. Kulkarni fe Frail I 
119981) . a nd through v ery lo ng baseline radio interfer- 
ometry ijTavlor et al. 1 12004(1 . At the sub-relativistic 
stage, the size was determi ned indirectly through mod- 
elling the radio spectrum llFrail. Waxman fc Kulkarni I 
[20001 iBereer. Kulkarni fc Frail Il2004j) . These observa- 
tions were used to determine the values of model param- 
eters that determine the flow pattern, E and n, under 
the assumption of high efficiency, / = 1. Modification, 
due to changes in the value of /, of the values of E and 
n are therefore allowed provided these changes do not 
modify the flow pattern. 

We demonstrate here that the velocity field v(f, t) as- 
sociated with the afterglow stage of the GRB explo- 
sion depends on the explosion energy E and surrounding 
medium density n only through the ratio E/n. If the 
surrounding medium density is not uniform, then for a 
given functional dependence g(r) of the density on coor- 
dinates, n(r) — nog(r) where hq is some normalization, 
then the velocity field v(f, t) depends on E and no only 
through the ratio E/uq. 

It is straightforward to demonstrate the validity of the 
above statements through examination of the hydrody- 
namic equations. These may be written as 

d^ v = 0, dp{nu^) = 0, (1) 

where u M is the four-velocity and T MI/ is the energy- 
momentum tensor. For the GRB afterglow flows, it is 
appropriate to assume an ideal fluid flow, i.e. T^ v — 
rf^p + (u^u u / c 2 )(p + e) where p and e are the fluid 
pressure and (proper) energy density respectively, and 
an ideal gas equation of state, e = nm p c 2 + (7 — l)~ l p 
where 7 is the adiabatic index (7 = 4/3,5/3 for rela- 
tivistic, non-relativistic particles respectively). It is now 
evident that if {u^(f, t),n(r, t),p(r,t)} is a solution of 
the flow equations, then multiplying the density and 
pressure by a constant K and leaving the velocity field 
unchanged provides another solution of the equations, 
{u M (r, t), n' — Kn(r, t),p' — Kp(r, t)}. Since the energy- 
momentum tensor of the new solution T'^ v is related to 
the energy-momentum tensor of the original solution T^ v 
by T'^ = KT^, the energy of the flow in the modified 
solution is larger by a factor K compared to the energy 
of the flow in the original solution. 

The argument given in the previous paragraph proves 
the statement that v(f,t) associated with the afterglow 
stage of the GRB explosion depends on the explosion 
energy E and surrounding medium density n (no) only 
through the ratio E/n {E/uq). Since this argument is, 
however, rather abstract, it may be useful to examine in 
some detail how the afterglow flow is affected at various 
stages as E and n are changed. This examination will 
also be useful for the discussion of § 12.21 For simplicity, 
we assume in the following discussion a uniform density 
of the surrounding medium (It is straight forward to gen- 
eralize the discussion to a non-uniform density). 

Let us first consider the flow associated with a spher- 
ical, relativistic blast wave. When the shock radius 
R is sufficiently large, compared to the size of the re- 
gion of initial energy deposition, the flow becomes self- 
similar, with shock lorentz factor determined by energy 
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conservation, E oc 7 2 ni? 3 im plying 7 oc {E/n) 1 ' 2 R 3 ' 2 
ijBlandford fc McKee 1 1197^1 . The self-similar flow is 
therefore completely determined by the ratio E/n. 

The early afterglow, on minute time scale, is pro- 
duced at the onset of the interaction of relativistic GRB 
plasma with the surrounding medium. At this stage, 
the highly relativistic plasma ejected by the GRB en- 
gine with Lorentz factor 7$, the "fireball," drives a for- 
ward shock into the surrounding medium, and a reverse 
shock is driven back into the fireball and decelerates it. 
Once the reverse shock crosses the fireball plasma shell, 
the flow approaches the self-similar behavior described 
abov e. This transition stage take place (e.g. iWaxmanl 
2003) at a radius which is the larger of (i) the radius at 
which the self-similar Lorentz factor 7 oc (£ , /n) 1 / 2 i? -3 / 2 
drops below 73, and (ii) the radius at which the thickness 
of the shocked plasma shell in the self-similar solution, 
-R/7 2 oc (n/E)R 4 , exceeds the thickness A of the plasma 
shell ejected by the GRB. The transition radius depends, 
therefore, on E and n only through the ratio E/n. 

If the fireball is jet like, rather than spherical, then flow 
is well described as a conical section of a spherical fire- 
ball as long as the jet opening angle is 9j > I/7. In this 
case, E should be understood as the "isotropic equiva- 
lent energy" , the energy that would have been carried by 
the blast wave had it been spherically symmetric. When 
the fireball decelerates to 7 < l/fl,, the jet i s frequently 
assumed to expand sideways ijRhoads i flQQi)'). The con- 
dition 7 ~ 1/(9 j implies that the radius Rj at which side- 
ways expansion begins is given by (n/E)Rj oc 9 2 . Rj 
depends, therefore, on E and n only through the ratio 
E/n. 

During the stage of sideways expansion, the jet 
does not significantly propagate radially. Finally, af- 
ter the stage of si deways expansion, the flow be- 
comes sub-relativistic ifFrail. Waxman fc Kulkarni 120001 
iLivio fc Waxmanl 12000(1 . and, if it becomes a spher- 
ical blast wave, the renewed expansion is described 
by the Sedov-von Neumann- Taylor solutions. At this 
stage the time dependence of the shock radius is R oc 
(E^/ny/H 2 / 5 (e.g. Chapter XII of Zel'dovich & Raizer 
2002). Here t is the time and the true energy of the explo- 
sion, corrected for the jet-like geometry, is Et = E9 2 /2. 
Thus, at this stage as well the flow depends on E and n 
only through the ratio E/n. 

2.2. Radiation 

Afterglow observations at all stages of flow evolution, 
from the non self-similar onset of fireball interaction 
with surrounding gas, through the self-similar expansion 
phase and subsequent jet expansion phase (if present), 
and including the final sub-relativistic phase, are consis- 
tent with synchrotron emission of radiation from elec- 
trons accelerated to a distribution of the type shown in 
fig- n with high efficiency, / si, Under the assumption 
/ = 1, the values of model parameters, {E, no, cb, e e } (as 
well as 9j and p), are determined. Let us now consider 
what modifications are introduced by allowing / <C 1. 
We argue that the emission of radiation from shock ac- 
celerated electrons from a flow with parameter choice 
{E' = E/f,n' = n //, 4 - fe B ,e' e = fe e ,f < 1} (and 
9j = 9j, p' = p) is similar to that obtained for the pa- 
rameter choice {E, no, es, e e } and / = 1, for any / in the 



range m e /m p < / < 1. 

Let us first consider the velocity fields of the two flows. 
The flow pattern v(r,t) in the modified, / < 1, flow 
is similar to that of the / = 1 flow, since the energy 
and density have both been increased by the same fac- 
tor 1//, leaving the ratio E/n unchanged (and since 
9 j = 9j). Next, we note that the magnetic field dis- 
tributions in the two flows are similar. As explained in 
the previous section, the energy density in the modified 
flow is larger than that in the original flow by a factor 
1//, e'(r,t) = e(r,t)/f. Decreasing the magnetic field 
equipartition fraction by a factor /, e' B = fes, ensures 
that the magnetic field energy density is similar in both 
flows. 

Finally, we argue that the density and energy distri- 
bution of accelerated electrons is the same in both flows. 
The number density of accelerated electrons is identi- 
cal in the two flows: The number density of electrons is 
larger in the modified flow by a factor 1// compared to 
that in the original flow, n' = no/ f, but only a fraction 
/ of the electrons in the modified flow are accelerated. 
The total energy density in electrons is also the same in 
both flows, since e' e e'(f,t) = e e e(r,t). The fact that the 
electron energy density and accelerated electron density 
is similar in the two flows does not ensure that the energy 
distributions of accelerated electrons are similar, since in 
the modified flow some part of the electron energy density 
is carried by the non shock-accelerated electrons. This 
part is small, however, and therefore the energy distri- 
butions of the accelerated electrons are similar in both 
flows, as long as m e /m p < f. 

To see this, we note that during the relativistic phase 
of expansion, the characteristic Lorentz factor -f e o of ac- 
celerated electrons, e e o = 7 e ow e c 2 , is approximately de- 
termined by the relation (1 — f)jm e c 2 + fj e om e c 2 — 
e' e -ym p c 2 — fe e ^m p c 2 . Since afterglow observations im- 
ply e e ~ 1 for / = 1, as long as m e /m p < f we have 7 e o ~ 
e e "fm p /m e independent of /. During the sub-relativistic 
regime, j e o is determined by the relation (1 — f)m e v 2 /2 + 
fleorn e c 2 = e' e m p v 2 /2 = fe e m p v 2 /2. Here too, j e0 is 
independent of /, approximately given by j e om e c 2 « 
e e m p v 2 /2, as long as m e /m p < /. It is important to note 
here, that since we have several examples where afterglow 
observations cover bo th relativistic and sub-relativistic 
evolution phases (e.g. iFrail. Waxman fc Kulkarni 1120001 
iBereer. Kulkarni fc Frail Il2004|) . the efficiency f should 
be similar at both stages. This independence of / on 
7 is not necessarily surprising, since / may be, e.g., a 
function oim e /m p alone. 

Since the two flows have similar velocity fields, similar 
magnetic field energy distributions and similar acceler- 
ated electron distributions, the afterglow radiation emit- 
ted by the accelerated electrons is similar for the two 
flows. The presence of a non shock-accelerated electron 
population may, however, modify the radiation pattern. 
This issue is discussed in the following section. 

3. SIGNATURES OF LOW EFFICIENCY 

Consider the possible presence of a large number of 
"thermal" electrons that entered the shock at energy 
7m e c 2 as measured in the shock frame. Assume that 
they are heated somewhat to a typical energy of n^m e c 2 , 
where n <C m p /m e . Here n is a parameter that expresses 
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our ignorance of the plasma physics that governs the 
electron heating beyond the energy 7m e c 2 , which the 
electrons bring into the shock from upstream. The pres- 
ence of a large population of these "thermal" electrons 
at energy jr]m e c 2 <C jrripC 2 (or at ~ rjm e v 2 <C m p v 2 ) 
may affect the emitted radiation by producing a new 
component of emission, or by producing a large syn- 
chrotron self-absorption optical depth, thus suppressing 
the emission from accelerated electrons. Since the en- 
ergy distribution of the non shock-accelerated electrons 
does not extend (by definition) to energies 3> r yrjm e c 2 
(or 3> nm e v 2 ), they may affect the emitted radiation 
at frequencies v < u m , where v m is the characteris- 
tic synchrotron emission frequency of electrons of en- 
ergy r]jm e c 2 (or rjm e v 2 ). This frequency is lower by 
a factor (r)m e /m p ) 2 than the characteristic synchrotron 
emission frequency v m of accelerated electrons at energy 
7m p c 2 (or m p v 2 , in which case the self-absorption is a 
line). Since the time dependence of the characteristic 
synchrotron emission frequency of accelerated electrons 
typically behaves as v m ~ 10 18 (t/100s)- 3 / 2 , i.e. peaks at 
X-rays on minute time scale and drops below the optical 
on a time scale of 10 hr, the characteristic synchrotron 
frequency of the non shock-accelerated electrons drops 
from ~ 300?7 2 GHz on minute time scale to ~ Q.irj 2 GHz 
on 3 hr time scale, 

v m w l?7 2 (t/lhr)- 3/2 GHz. (2) 
This implies that the existence of a large population of 
non shock-accelerated electrons can be constrained only 
through early radio observations, at < r/ 4 ^ 3 hr delay. As- 
suming a spread of several times the thermal energy, so 
that 77 > 3, the emission and/or absorption of the ther- 
mal electrons could be caught with radio follow-up obser- 
vations within several hours of the GRB. With sufficient 
preparation, radio follow up observations may be carried 
out on minutes time scale (D. Frail, G. Taylor, private 
communication) . 

The specific emissivity j m of the non shock-accelerated 
electrons at v m is larger than that of the accelerated 
electrons at v m by a factor 1// (since their number is 
larger by this factor). Thus, if the optical depth at 
v m is small, the radio intensity produced by the non 
shock-accelerated electrons would be larger by a fac- 
tor l/f 1 than the ~ 1 mJy peak intensity char- 
acteristic of the accelerated electrons for cosmological 
GRBs. The synchrotron self-absorption optical depth 
at the peak frequency can be estimated using Kirchoff 's 
law, f m oc j m /v^T where the effective temperature is 
T = njm e c 2 . From this relation, we find that the ra- 
tio of f m to the optical depth r m at v m is (for small 
/) T m /T rn rs (mp/r)m e ) f ■ For the population of 
accelerated electrons, r m = (v a /v m ) 5 ^ 3 where the self- 
absorption frequency v a ~ 1 GHz and inde pendent of 
time fo r expansion into uniform medium (e.g. iWaxmanl 
1997b). Combining these relations we have, for expan- 
sion into uniform medium and small /, 

/8 f? - 6 /- 1 no(t/lhr) 5 / a . (3) 



We have kept here the dependence on the ambient 
medium number density, n = 10 n riocm -3 , mainly in or- 
der to allow a simple generalization to the case of ex- 
pansion into a non uniform medium. For expansion 



into a wind, n oc t 1 , and for typical wind parameters 
no » l(t/l day)" 1 (|Livio fc Waxman 11200(1 . All the re- 
sults given here can thus be applied to the wind case by 
using uq — l(t/l day)" 1 (note, that eq. (J2J is valid for 
any density, i.e. has no dependence on n). 

The optical depth at f m is larger than unity for t > t a , 
where 

t Q «10-V 7 V/ 2/5 hr. (4) 

At t < t ai the self-absorption frequency v a , where 
the optical depth due to the thermal electrons is unity 
{f m {i> m /v a ) 5 / 3 = 1), is 



t>p 3/5 



Ly-3/5 GHz 



(5) 



Finally, the specific intensity at i> n 

fm ~ f^fm x min [1, l/f m ] 
1/2. 



■f 



-1 "0 



frr 



lmjy 



is given by 



l,(t/t a )-V 2 mJy. (6) 



Here, f m ~ 1 mJy is the peak intensity characteristic of 
the accelerated electrons for cosmological GRBs. 

The presence of a significant number of non shock- 
accelerated electrons (/ <C 1) is therefore expected to 
lead to a large self-absorption optical depth at frequen- 
cies v < P m , strongly suppressing the radio flux at these 
frequencies at early time. A sharp rise in the flux at fre- 
quency v is expected at t m 77 4 / 3 (^/lGHz)~ 2//3 hr, as v rn 
drops below v. The ratio of fluxes obtained in the pres- 
ence and in the absence of a non-accelerated electron 
population, f v j '/„, is given for frequencies lower than 
the self-absorption frequency of the non-shock acceler- 
ated electrons, v < min(f' m , v a ), by the following argu- 
ment. At frequencies where the optical depth is large, 
the intensity is proportional to v 2 T where T is the ef- 
fective electron temperature. For the non accelerated 
electrons, T = r]"fm e c , and, in the absence of electron 
cooling, T = 7 e o?n e c 2 , so that T/T w r)m e /m p . At early 
times, the cooling time of accelerated electrons is short 
compared to the dynamical (expansion) time, and these 
electrons lose energy and accumulate at lower Lorentz 
factor, 7c, where the cooling time is comparable to the 
dynamical time. At this energy these electrons radiate 
synchrotron photons at frequency v c — {^cheo) 2 v m , and 
for typical model parameters 



— w 10n (t/lhr)- 



(7) 



Thus, for v < min(£' f , 
given by 

fu _ yrn e 
U ~ 



, Vq) the flux suppression factor is 



■ max 



1, 



1/2' 



1,1- 



5/3' 



(8) 



The last term on the rhs accounts for the fact that 
f v oc v x l z (rather than f v oc v 2 ) for v > v a . Cooling of 
electrons increases v a by a factor "f e o/"fc — (y m j Ve) 1 / 2 , 
compared to the case where cooling is unimportant, 
which implies 



3/5 
n max 



l,3ng /2 (i/lhr) 



-1/2 



GHz. 



(9) 



The presences of a large number of non shock- 
accelerated electrons may be detected through their ra- 
dio emission only if this emission takes place in an op- 
tically thin regime, i.e. at frequencies i> a < v < v m 
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for v a < v m . Examining eqs. JSJ and (J2J), we find that 
v a < v m is possible only for r\ 3> 1 and moderate / . 
For rj = 10 and / = 10 _1 , for example, the flux at 
v m RJ 400(n/10) 2 (i/0.4hrT 3 / 2 GHz is - lO^O/)" 1 mJy 
up to t rj 0.4(n/10) 2 (10/) 2 / 5 hr. For v a < v < v m and 
v > v a we have 



./V 



-2/3 



-1/3' 



(10) 



Finally, it is useful to give an estimate of the amplitude 
of the radio flux typically expected on the relevant time 
scale. For v > v a , the flux expected (in the absence of 
non shock accelerated electrons) is 



-1/3' 



1/3 



fn 



:30n ( 



i/6 fn 



x max 



lmjy 

v r 



(k^e)" 3 "/ 11 "'"** 



-1/3' 



(11) 



4. DISCUSSION 



We have shown that current afterglow observations do 
not allow one to determine the efficiency of electron ac- 
celeration in GRB shocks, i.e. to determine the fraction / 
of electrons that are " injected" to participate in the pro- 
cess of shock acceleration. While afterglow observations 
imply that some fraction, /, of the electrons' population 
is accelerated to a characteristic energy e e Q comparable 
to the post shock proton temperature, e e o RJ jrripC 2 for 
relativistic shocks of Lorentz factor 7 or e e o RJ m p v 2 /2 
for non-relativistic shocks of velocity v, a large fraction, 
1 — / ~ 1, of the electron population may be "left be- 
hind" at low energy comparable to the kinetic energy of 
the electrons propagating to the shock, 'ym e c 2 for rela- 
tivistic shocks or m e v 2 /2 for non-relativistic shocks. The 
resulting electron energy distribution is qualitatively de- 
scribed in figure n Currently testable afterglow predic- 
tions of a model with parameter choice {E' = E/f,n' = 
no/ f, e' B — feB,t' e — fte, f < 1} are similar to those ob- 
tained for the parameter choice {E, no, es, e e } and / = 1, 
for any / in the range m e /m p < / < 1. This implies 
an uncertainty of factor m e /m p in the determination of 
model parameters. Afterglow observations do not con- 
strain, for example, the values of E and e e , but rather 
the values of fE and e e /f. Note, that the value of e e o is 
independent of / (and equals rj jm p c 2 or rj m p v 2 /2). 

The existence of non shock-accelerated electrons will 
strongly affect the predicted radio emission on short, 
< in 4 / 3 nr; time scale. Here, r]~/m e c 2 (or r]m e v 2 /2) is 
the characteristic energy of non shock-accelerated elec- 
trons (77 <C m p /m e ). For / « 1, a large self- absorption 
optical depth at v m rj l?y 2 (£/lhr)" 3 / 2 GHz [eq. would 
lead to strong suppression of the radio flux at lower fre- 
quencies [eq. ©]. As v m drops below an observed fre- 
quency v : at t rj n 4 / 3 (^/lGHz)~ 2 / 3 hr, the optical depth 
at this frequency drops below unity, and a sharp bright- 
ening is expected. For i] ^> I and moderate / , the 
existence of a large population of non shock-accelerated 
electrons may be identified through their radio emission 



[see eq. ijTUj l and the discussion preceding it]. For rj 3> 1, 
the modification of radio emission due to the presence 
of non shock-accelerated electrons will persist over time 
scales significantly larger than 1 hr. The radio signature 
of these thermal electrons could test for their presence 
at levels that are energetically insignificant by a large 
margin (even / _1 < I), and therefore otherwise incon- 
spicuous. 

It should be pointed out that afterglow observa- 
tions already provide interesting constraints on the ef- 
ficiency of electron acceleration. First, they require simi- 
lar efficiency / for both relativistic and sub-relativistic 
shocks, since several examples exist of afterglow ob- 
servations covering b oth relativistic and sub-relativistic 
evolution phases (e.g. iFrail. Waxman fc Kulkarni 1120001 
iBerger. Kulkarni fc Frail 1120 04). This independence of / 
on 7 is not necessarily surprising, since / may be, e.g., 
a function of m e /m p alone. Second, afterglow observa- 
tions imply that the energy of accelerated electrons is 
increased to a characteristic energy similar to the pro- 
tons post-shock temperature (with power-law extension 
to high energies). Finally, the value of / is limited to 
/ > m e /m p . Early, < I hr, radio observations will pro- 
vide more stringent constraints on the efficiency / (and 
will hence remove the degeneracy in determining GRB 
model parameters) . As mentioned in § |2 radio spectra 
can also be used to constrain the energy distribution of 
accelerated electrons at ene rgies below the ch aracteristic 
acceleration energy, jrripC 2 ( Wa xman"lll997b() . providing 
further constraints on the acceleration process. 

The total, beaming corrected, energy released 
in cosmological long duration GRB explosions 
is typically inferre d, assuming / = 1, to be 
10 51 5 erg (IFrail Waxman fc Kulkarnil 



E- 



T 



Bcrgcr, Kulka rni fc Frail 



2001: 



| Frail et al.T 

2004), with 



2000; 



2001; 



1200.11 1200-41) . with a spread 
in estimated values of roughly one order of magni- 
tude. Since afterglow observations do not constrain 
Et, but rather Er/f, the true explosion energies are 
E T ~ / _1 10 515 erg. For / < 1, explosion energies 
3> 10 erg would naively be inferred for many 
GRBs. The association of (at lea s t some) GRBs with 
supernovae iGalama et al. 1119984 iStanek et al~l 12003 
iHiorth et al. I l200.1t iBloom I 120031) suggests that the 
total energy is probably not much more than 10 515 erg, 
ruling out values of / <U 1. Using this argument to infer 
a conservative lower limit on /, the uncertainties in 
determining Et / f from afterglow observations should 
be considered. Uncertainties in determining E may 
arise from uncertainties in the determination of the 
observables {v m , F m ,v coo i,v sa ). The main uncertainty 
here is due to uncertainty in the determination of the 
self-absorption frequency, which is not known in many 
cases and determined in the best cases to within a 
factor of ~ 2, leading to uncertainty in E of a similar 

magnitude (since E oc Vsa^ ', e.g. iWiiers fc Galamal 

1999, Note, that the uncertainty in n oc i/j^ 6 and 

€b oc v sa is much larger). Moreover, it should be 
realized that afterglow models are highly idealized 
(e.g. assuming simple geometry) and various effects 
which are not taken into account (e.g. acceleration of 
pre-shocked plasma by a cosmic-ray precursor) may lead 
to systematic errors in estimates of model parameters. 
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Thus, estimates of the energy should be considered 
as order of magnitude estimates. Finally, the energy 
provided by the supernova to the GRB jet could be 
higher than that provided to the supernovae ejecta, 
whose energy is limited by neutrino cooling. Altogether, 
although / ~ 1 is suggested by the energy derived from 
afterglow observations, / > 1/30 should be considered a 
plausible conservative lower limit. 

Is there an a priori reason to suspect that / should 
be small? In the case of plerionic nebulae, such as the 
Crab, typical shock-accelerated spectra (p ~ 2.2) oc- 
cur above Lorentz factors of order e ~ 10 4 , and be- 
low that the spectra are much flatter, 1.3 > p > 1.1 
ijWeiler fc Panagia Il978j) . Curiously, the low energy end 
of these spectra goes well below 7m e c 2 (where the bulk 
Lorentz factor of the pre-shock wind can be estimated 
knowing the total number of electrons and the total en- 
ergy that have been deposited into the nebulae by the 
wind). These low energy electrons do not increase the 
total energy requirements, but clearly comprise most of 
the electrons by number. This raises the question of how 
most of the electrons in the nebula can have less en- 
ergy than they had flowing into the shock, and strongly 
suggests some sort of shock mediation mechanism that 
redistributes their energy in the form of a hard power 
law. 

Suppose the same sort of low energy spectra were ob- 
tained below 777i p c 2 in GRB blast waves. We can ex- 
press these low energy electron populations as / _1 (?7) — 
(rjm e /m p )~ p+1 . Basically / _1 (?7) is the number of low 
energy electrons at energy rjm e relative to those at m p . 
If plerion-like low energy spectra were to obtain in GRB 
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post-afterglow-shock plasmas, it would give values of / 
of more than 1/30 at rj > 1 but would violate the con- 
straint d\nn/d\nE > that seems to exist at least 
for some GRB afterglows. These current limits ex- 
tend from 1/30 < rjm e /m p < 1. Exploring the region 
r}m e /m p < 1/30 at I^ghz GHz will, by equation (2), re- 
quire radio follow up observations within (10 3 ' 5 /fGHz) 2//3 
hr. The signature of a thermal population of electrons 
would be a) a radio "blackout" due to the low bright- 
ness temperature of the thermal electrons, followed by 
b) a pre-brightcning, as the emitting area increases, fol- 
lowed by c) a steep decline, as the emitting frequency 
of the thermal electrons passes below the observed fre- 
quency. All these stages precede the expected rise as- 
sociated with the eventual passage of v m through the 
observing frequency. 

In summary, early observations of emission and self- 
absorption in GRB afterglows can provide a diagnostic of 
the low energy electron spectra in GRB afterglow shocks. 
This might either change our understanding of them or, 
at least, close some important loopholes in afterglow the- 
ory. 
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